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ABSTRACT 

We investigate the star formation history of the central regions of four Luminous Compact 
Blue Galaxies (LCBGs) at intermediate redshift using evolutionary population synthesis tech- 
niques. LCBGs are blue {B — V < 0.6), compact (/is < 21.0 mag arcsec"^) galaxies with absolute 
magnitudes brighter than -17.5. The LCBGs analyzed here are located at 0.436 < z < 0.525. 
They are among the most luminous {Mb < —20.5), blue {B — V< 0.4) and high surface bright- 
ness {fiB < 19.0 mag arcsec"^) of this population. The observational data used were obtained 
with the HST/STIS spectrograph, the HST/WF/PC-2 camera and the HST/NICMOS first cam- 
era. We have disentangled the stellar generations found in the central regions of the observed 
targets using a very simple model. This is one of the first times this is done for compact galaxies 
at this redshift using Hubble Space Telescope data, and it provides a comparison bench for the 
future work on this kind of galaxies using instruments with adaptive optics in 10-m class tele- 
scopes. We find evidence for multiple stellar populations. One of them is identified as the ionizing 
population, and the other one corresponds to the underlying stellar generation. The estimated 
masses of the inferred stellar populations are compatible with the dynamical masses, which are 
typically 2-lOxlO^M0. Our models also indicate that the first episodes of star formation the 
presented LCBGs underwent happened between 5 and 7 Gyr ago. We compare the stellar popu- 
lations found in LCBGs with the stellar populations present in bright, local H II galaxies, nearby 
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spheroidal systems and Blue Compact Dwarf Galaxies. It turns out that the underlying stellar 
populations of LCBGs are similar yet bluer to those of local HII galaxies. It is also the case 
that the passive color evolution of the LCBGs could convert them into local Spheroidal galaxies 
if no further episode of star formation takes place. Our results help to impose constraints on 
evolutionary scenarios for the population of LCBGs found commonly at intermediate redshifts. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift — galax- 
ies: stellar content 



Introduction 



Luminous Blue Compact Galaxies (LCBGs) are luminous {Mb brighter than -17.5), blue {B — V< 0.6), 
and compact (/is < 21.0 mag arcsec"^) systems experiencing an intense episode of star formation, with 
typical star formation rates ranging from 1 to 5 MQyr~^. Although these star formation rates are common 
in other starburst systems, the star formation episode takes place in a very small volume, allowing these 
compact objects to be seen from cosmological distances. LCBGs are believed to contribute up to ~50% of 
the star formation rate density in the universe at z ~ 1 (Guzman et al. 1997). Their rapid evolution from 



those times till the present day makes them responsible for the observed evolution of the luminosity function 
( Lilly et al. ||1995 ) as well as the evolution of the star formation rate in the universe. 



Despite their importance, their relation to local galaxies is poorly known. The works presented in 



Guzman et al. (19971; Phillips et al. (1997); Hoyos et al. (2004) give a glimpse of the similarities between 



LCBGs and the local HII and Starburst Nuclei populations. HII galaxies are a subset of Blue Compact 
Galaxies in which the spectrum is completely dominated by a component, present almost everywhere within 
the galaxy, resembling the emission of a HII region (Sargent & Searle 1970). Blue Compact Galaxies 



(BCGs) were first identified by Zwicky ( 1965 1 as faint star-like field galaxies on Palomar Sky Survey plates. 
BCGs often present a emission line spectrum and a UV excess. Starburst Nuclei were introduced in |Balzano | 



( 1983 ); they show high extinction values, with very low [Nil] A 6584/Hq; ratios and faint [OIIIJA5007 emission. 
Their Ha luminosities are always greater than 10^ solar luminosities. One of the main ingredients in any 
theory attempting to explain the final destiny of LCBGs is bound to be an accurate knowledge of the stellar 
generations already existing in these distant sources. The similarities found between LCBGs and many local 
HII galaxies make it possible that the star formation histories of both types of objects are very similar. 
Even though H II galaxies were originally thought to be truly young objects, experiencing their very first 
star formation episodes ( Searle fc Sargent|1972 Sargent & Searle 1970), this hypothesis is no longer tenable 

1997, ,Telles, Melnick fc Terlevich |p97| . 



after very deep CCD imaging of such objects ( Telles & Terlevich 



or other very similar sources like the Blue Compact Dwarf Galaxies (BCD). BCDs are similar to BCGs or 
HII galaxies, but are an order of magnitude less luminous. The underlying stellar populations of BCDGs 
are also older than those of BCGs. ( Papaderos et al.| 1996a Loose fc Thuan|1986 or Papaderos et al.|1996b l 
which show clear evidence of an older extended population of stars. The properties of this more extended 



component has been further studied photometrically (see, e.g. Caon et aL]|2005 Telles, Melnick & Terlevich 



1997). The current view on local H II galaxies is that they are a mixture of three populations. The first 
of them is the youngest population, responsible for the observed emission lines. The second population is 



^ Based on observations obtained with the NASA/ESA Hubble Space Telescope through the Space Telescope Science Institute 
which is operated by Association of Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS 5-26555 
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an intermediate age generation, older than 50Myr but younger than 1 Gyr. This population is unable to 
ionize the H I gas, and its luminosity is increasingly dominated by the contribution of giant stars, and by 
the contribution of some Asymptotic Giant Branch stars. The third population is the oldest one, older than 
IGyr. This population generally dominates the stellar mass by a large factor, and it also occupies a larger 
volume. 

The main target of this study is to characterize the stellar populations present in intermediate redshift 
LCBGs by combining the spectral and color information available from these data with evolutionary popu- 
lation synthesis techniques. The modelled stellar populations are then compared to the stellar populations 
found in HII galaxies and the smaller yet similar Blue Compact Dwarf galaxies. It is also interesting to 
study the future evolution of the predicted stellar populations, to see whether or not they can become dwarf 
elliptical systems (dE, or Sph), as suggested by Koo et al. (19951; Guzman et al (19981 or, alternatively, 



they can be the bulges of spiral galaxies still in their formation process ( Hammer et al. 
fc Zaritsky ||1999[ [Guzman e^"al]|1998| [Phillips et aT]|1997[ ) . 



2001 Kobulnicky 



The observations and analysis techniques together with the different results they yield are described 
in section |2] Section |3] presents the data analysis. The models that have been constructed to explain the 
observations are explained in section [4] The LCBGs presented in this work are compared to relevant samples 
of local galaxies in section [5] The summary is given in section [6] 

The cosmology assumed here is a flat universe with — 0.7 and ^matter — 0.3. The resulting 
cosmological parameters are: Ho=70 km s^^Mpc^^, qo=-0.55. Given these values, 1" corresponds to l.Skpc 
at z = 0.1 and to 5.8 kpc dX z — 0.45. 



2. Observations and Data Reduction. 

In this paper we present long-slit spectra, deep optical V (F606W) and I (F814W) images, and infrared 
F160W images of four LCBGs with redshifts between 0.436 and 0.525. The data were obtained using the 
STIS instrument, the WF/PC-2 camera and the NICMOS first camera. Target names, redshifts, look-back 
times, B absolute magnitudes, exposure times, half light-radii, velocity dispersions and virial masses within 
i?e, as well as an identifying numeral for the observed objects can be seen in table [l] 

The STIS long slit spectra are imaged on a SITe 1024x1024 CCD, with 0.05" square pixels operating from 
~2000 to 11000 A. Several instrumental configurations were used. The instrumental setup used to observe 
objects 1, 3 & 4 used the low resolution grating G750L. This instrumental setup provides a dispersion of 
4.92 A per pixel or about 190 km s ^ per pixel which, in combination with a slit 0.5" wide, gives a spectral 
resolution of about 29A(FWHM). The spectral coverage provided by this setup is about 5000 A. The central 
wavelength is 7750 A. The higher resolution grating G750M was also used for objects 3 & 4. It provides 
a dispersion of 0.53 A per pixel or about 35 km s~^ per pixel. The slit width used was 0.2" wide, thus 
giving a spectral resolution of about 2.lA(FWHM). The spectral coverage provided by this setup is about 
500 A. These spectra are used to study the gas-phase velocity field, and the results from this work are 



presented in Bershady et al. (2005). Object 2 was observed using the lower-resolution grating G430M and 
the high-resolution grating G750M. The G430M grating has a dispersion of 2.7 A per pixel and the slit 
width used is 0.2". The central wavelength is 4300 A, and the spectral coverage is 2800 A. The data were 
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Table 1. Log of HST/STIS, HST/WF/PC-2 and HST/NICl observations. 



Parameter. 


bA57-7042. 


bA57-5482 


SA57-10601. 


Hl-13088. 




13 07 26 3 


13 09 08 8 


13 08 47 8 


17 20 19 67 


DEC^ 


29 18 25 7 


29 1 5 57 7 


29 23 41 1 


50 01 04 7 


Obs Date. 


2002-05-05 


2002-04-09 


2000-07-05 


2000-11-06 


G750L Exp. Time. 


9898 




4696 


4932 


G750M Exp. Time. 


2760 


18106 


5571 


5947 


G430L Exp. Time. 




2790 






G430M Exp. Time. 




7709 






WF/PC-2 F606W Exp. Time. 


700 


700 


700 


700 


WF/PC-2 F814W Exp. Time. 


900 


900 


900 


1000 


NICMOS F160W Exp. Time. 


2560 


1280 


5632 


1792 


Obj ID. 


1 


2 


3 


4 


z. 


0.525 


0.453 


0.438 


0.436 


tlook-back{Gyr) 


5.2 


4.7 


4.6 


4.5 


Mb 


-20.6 


-20.8 


-20.4 


-21.1 


Re'' 


1.4 


1.5 


1.9 


2.4 


IBB 


22.3 


21.6 


21.6 


21.0 


r 

'-e 


0.25 


0.29 


0.36 


0.47 


Vel. Disp. (km s^"'^ 


120±10 


60±6 


42±4 


47±4 


logM/Mo 


10.7 


10.1 


9.7 


9.8 



^Coordinates given for equinox J2000. 
'^Half-light radii in kpc. 
■^Half-light radii in arcseconds. 



''Dynamical masses within i?e calculated as in Hoyos et al. (20041. Uncertainties are 
20%. 



Note. — All quantities given were calculated assuming go = —.55 and Hq — 70km 
s-'Mpc-'. 
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reduced using the STSDA^ package within IRAF ^ For a complete account of the reduction procedures, 



uncertainties and detailed descriptions of the spectra, see iHoyos et al. (2004). It is worth mentioning that 



STIS spectroscopic observations using the G750L and G750M gratings are affected by fringing in the red 
end of the spectrum. This makes it neccessary to use specifically designed flats (called "contemporaneous 
flats" ) in order to circumvent this issue. Thes use of this technique, which is implemented as a task within 
STSDAS raises the S/N by around 15%. This is crucial for these observations. The proposal IDs for these 
observations are G08678 and G09126. 

The images used were taken with the WF/PC-2 instrument, with the LCBGs centered in the Planetary 
Camera, a 800 x 800 CCD, with 0.046" square pixels. Two broadband filters were used, the F814W filter 
and the F606W filter. For each object and filter, two ~ 400s exposures were obtained, to allow for cosmic 



ray rejection. The details of the reduction and analysis procedures can be found in Guzman et al (19981. 
The proposal ID for these data is G05994. 

The infrared frames were obtained with the NICl camera of the NICMOS instrumeni|f] This camera 
uses a low- noise, high QE, 256x256 pixel HgCdTe array, with 0.0431" square pixels. The filter used was the 
F160W, which has a bandpass of around 4000A. For each object, several different exposures were collected 
and combined to produce a final mosaic. The proposal ID is G07875. 

The STIS data can only probe the central regions of these galaxies because the outermost regions of 
these sources are too faint to be studied spectroscopically with HST. Specifically, the STIS data can be used 
to study the line-emitting region, where the starburst is taking place and most of the gas is ionized. The 
line-emitting region of the observed galaxies has a radius between 1.4 and 2.0 kpc. This is a little bit less 
than one effective radius. On the other hand, the WF/PC-2 images and NICMOS frames are deeper, and 
the size of the observed galaxies in the images is much larger than the size of the line emittiong region in 
the STIS spectra. The analysis presented in this paper focuses on this particular central region where the 
STIS spectra spatially overlaps with the WF/PC-2 images and NICMOS frames. This region has a radius 
of 1.7 kpc, on average. 



3. Results and Analysis. 

3.1. Emission-Line and Continuum Distribution. 

It is interesting to investigate whether or not the continuum distribution and the emission-line dis- 
tribution differ in the STIS spectra. This will help determine where the starburst has taken place. This 
information can be helpful to understand the starburst triggering mechanism, and to construct models for 



these objects. Here, the analysis presented in Hoyos et al. (2004) is followed, although it is presented in 



a much simpler way. Furthermore, only objects 1 & 2 are presented here since the other two sources were 



dealt with in Hoyos et al. (2004). The procedure is briefly reviewed. 



Using the IRAF task fitl d, a "continuum frame" was made by fitting 3 cubic spline pieces to every line 



^ (Space Telescope Science Data Analysis System) from the Space Telescope Science Institute, operated for NASA by AURA 

^ Image Reduction and Analysis Facility. Distributed by the National Optical Astronomy Observatories, which is operated 
by AURA (Association of Universities for Research in Astronomy, Inc) under cooperative agreement with the National Science 
Foundation. 



^Before the cryo cooler was installed. 



-6 - 



in the spectra, and then, a "Hne frame" was constructed by subtracting the "continuum frame" from the 
original 2-D spectroscopic frame. The spatial distribution of the continuum is obtained by column-averaging 
the whole "continuuum frame," and the spatial distribution of the emission lines are obtained by column 
averaging the "line frame" over the columns that span 1.25x the FWHM of the line. The spatial profiles for 
the two new objects studied are given in figure [l] 

Figure [T| shows that: 

1. The spatial distribution of the ionized gas and underlying stellar population in object 1 is very similar 



to the distributions observed in objects 3 and 4 (Hoyos et al. 20041: The nebular component light 



centroid does not coincide with the continuum centroid although the peaks of the different components 
clearly line up. The two peaks observed at +1.5 kpc and -2.7 kpc in the H/3 spatial profile are not 
real. This object is very similar to the cometary BCDGs. In this particular case, the spread of the 
continuum is about twice the size of the line-emitting region. In the optical images, this object shows 
a conspicuous tail, and a very bright knot. It is also seen that the spatial extent of the two emission 
lines used to trace the ionized gas is very similar. 

2. For object 2, the diferences in the centroid positions for the continuum and the line-emitting region are 
very small when compared to the spread of both. In this case, the ionized gas can be said to lie in the 
center of the optical galaxy. This is to say that the line emitting region coincides with the continuum 
light centroid. The true barycenter of these galaxies is not known. It is also interesting to note that, 
in this particular case, the distributions of each component are very noisy even for the continuum, and 
that the distributions of H/3 and [O III]A5007 are very similar to each other. 



As in Hoyos et al. (20041, the H/3 and [O III]A5007 lines were chosen to probe the ionized gas phase 
because [O III]A5007 and H/3 are very close to each other in wavelength space and therefore they are affected 
by a similar extinction. For this reason, differences in their spreads can then reveal real differences in their 



distributions and not just a reddening effect. As it was pointed out in Hoyos et al. (20041, the observed 
differences in the spread of Ha and H/3 or [O III]A5007 and [O II]A3727 can be accounted for by a moderate 
extinction, within the measurement errors, making it impossible to detect extinction or excitation gradients 
unambiguously. Finally, the inclusion of these two new objects to the sample of LCBGs with high quality 
spatially resolved spectroscopy allows us to give a first interval of the frequency of off-centei[^ starbursts in 
LCBGs. Since it results that three out of six LCBGs with STIS spectroscopy happen to have their starburst 
displaced from their centers, all that can be said is that the fraction of LCBGs with shifted starbursts is 
50l;2o9o- More observational work is needed to asses the frequency and impact of this issue in the evolution 
of LCBGs. 



3.2. Line Ratios, Reddening, Equivalent Widths and Metallicities. 

The observed flux of each emission line and the equivalent widths have been measured via gaussian fits 
using the STSDAS task ngaussfit, according to the method outlined in Hoyos et al. (2004). This task can 
fit to a set of x-y values the sum of a straight baseline, and one or several gaussian profiles. This task also 
gives approximate errors for the fitted parameters, calculated by resampling. The emission lines were fitted 
rather than summed because of residuals of bad pixels. 



^Here, the word "center" means the continuum light centroid. 
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Two different positions are considered for each galaxy corresponding to (a) the inner and (b) the outer 
regions of each object. Despite these names, both the (a) and (b) regions defined here are within the hne- 
emitting region. These two one-dimensional extractions are defined as to encompass approximately the same 
line luminosity. It is in these two regions that the following stellar population analysis is carried out. These 
apertures are near one eff'ective radius for the line-emitting region, as detected through the spectroscopic 
aperture. The linear sizes of both extractions are given in the second row of Table [2] in kpc. For position (a), 
both the FWHM and central wavelength are left as free parameters. However, for position (b), the FWHM 
is fixed to the FWHM derived for position (a) to decrease the number of free parameters and hence increase 
the goodness of the fit in poor S/N conditions. The [O II]A3727, H/? and both [O III] lines can be observed 
in all cases. However, athough Ha is detected for object 4, it is affected by many bad-pixels that render the 
line unusable. In the case of object 1, H7 is detected, too. 

The amount of extinction is here parametrized using the logarithmic extinction coefficient c(H/3). It 



was derived using standard nebular analysis techniques (Osterbrock 19891 using the Ha/H/3 Balmer line 



ratios. For object 1, the H7 line was used to derive the extinction coefficient, which agreed well with that 
derived from Ha and H/3 but had a smaller uncertainty. Since the observed galaxies are at high Galactic 



latitude, suffering galactic extinctions Ab < 0.11 mag which corresponds to a c(H/9)< 0.05 (Burstein & Heiles 



1982 Sandage & Tamman 19741, it was assumed that all the extinction occurs within the target objects. 



This approximation is further validated by the fact that the incoming light is redshifted and therefore the 
extinction will be much lower. Consequently the spectra have been dereddened in the rest frame of the targets 
according to a standard reddening curve dWhitford |p58| , assuming case B recombination theorethical line 
ratios. The use of different extinction curves would not significantly change the results since they do not 
differ from the one used here in the wavelengths of interest. The value of c(H/3) is determined independently 
for the inner and outer regions for objects all objects but Hl-13088, for which Ha is unfortunately not well 
measured. A global value of c(H/3) = 0.1 ± 0.05, typical for local H II galaxies (Hoyos & Di'az 20061 is 



adopted. No underlying stellar hydrogen absortion is detected, as the continuum is very low. 

Line equivalent widths were measured in the STIS spectra. It is possible to measure the line equivalent 
widths despite the fact that these spectra are emission-line dominated, if the error estimates are carefully 
estimated. This is explained in detail in [Hoyos et al. ( 2004[ |, where it is also possible to see four examples of 
the extracted spectra (see again [Hoyos et al.|[2004l Fig 2.). 

The measured extinction coefficients and de-reddened line ratios are given in Table [2] together with the 
measured rest-frame equivalent widths. It is seen that the extinction values are always very small. It can 
also be seen that the ionization ratio log [O ni]/[0 II] and both the log [O III]/[H II] and logi?23 numbers are 
always smaller in the outer zones of the observed galaxies than in the inner part of the line-emitting region. 
R23 is defined as the reddening-corrected ([O II]A3727-|-[0 III]A, A4959, 5007)/H/3 ratio. In some cases, the 
differences are larger than the derived errors revealing the existence of a gradient. For some other cases the 
differences are more dubious. However, even though it is not possible to claim that these gradients have 
been detected for all the individual objects studied, the fact that these differences between the (a) and (b) 
zones hold for the majority of cases allows to think that these gradients are real at least for typical LCBGs. 



3.3. The Oxygen Content. Metallicity Estimates. 



One of the main goals of any nebular analysis is to estimate the metallicity of the emitting nebula. 
There are many reasons to determine the metallicity of the ionized gas in a galaxy, but chief among them 
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is that it yields an up-to-date determination of the chemical composition of the most recent generation of 
stars born in the galaxy studied. 

The emission line spectrum from a star forming region depends strongly on the electron temperature 
and metallicity of the nebula. The ratio of the auroral line [O III]A4363 to the lower excitation lines 
[O III]AA4959, 5007 gives a direct determination of the electron temperature where 0+ and 0++ are the 
dominant species, but beyond certain metallicities, the increasing cooling leaves no energy to collisionally 
excite the upper levels. When this is the case, one must resort to determine the oxygen abundance empirically. 

In principle, the oxygen content can be determined from the intensities of the optical emission lines 
using the i?23 number, defined as the reddening-corrected ratio ([O II]A3727-|-[0 III] A, A4959, 5007)/H/3. This 
method was originally developed by Pagel et al. ||1979 (see also Pagel, Edmunds fc Smith ||1980[ [Edmunds 



fc Pagel]|1984 McGaugIi||1991 Kewley fc Dopitar||2002 among many others). The main problem with this 
empirical method is that the calibration is double- valued, and some a priori knowledge of the metallicity 
range is needed in order to solve this degeneracy. Since the S/N of the spectra used is not very large, no 
other lines that could help break this degeneracy are detected. Thus, the lack of this initial knowledge about 
the oxygen content remains, and the metal content of the galaxies studied cannot be cleanly determined. 
Instead, the metallicity values that result from the use of the Pilyugin, L. S. ( [2000 ') calibration for both 
the lower and upper branches are calculated. The metallicity estimates calculated using the Pilyugin, L. ST| 
I ( 2000 1 formulae are given in Table |2] Unfortunately, the measured values for log i?23 squarely place our 
objects in the turnaround region of the 12 -I- log (0/H)-log i?23 plot. In this regime, the only thing that 
can be safely said is that the metallicity is between 7.7 and 8.4. Therefore, the analysis can only safely 
conclude that these LCBGs have an oxygen content greater than 7.7, but lower than the solar value of 8.69 



( AUende-Prieto, Lambert & Asplund 2001) It is also worth mentioning that, with the available data, it is 



not possible to unambiguosly detect any gradient in the metallicity of these objects since the derived values 
of i?23 imply that the oxygen content of both zones can lie in a very wide rang^ It can also be concluded 
that the i?23 method is not best suited to determine the gas-phase metallicity of LCBGs because the i?23 
diagnostic is likely to fall in this regime. 

However, one initial metallicity estimate can be provided by the ratio [N II]A6584/Hq! as suggested by 



Denicolo, Terlevich & Terlevich (2002 1. Since [N II]A6584 can not be measured in the spectra analyzed here. 



as it is below the detection limit due to poor S/N, it is not possible to use the spectra used here directly for 



this task. The only available information are the Guzman et al. (1996) W.M. Keck Telescope [N II]A6584/Ha 
measurements for a set of very similar sources from the same parent sample. These measurements consist of 
very high quality spectra obtained using the High REsolution Spectrograph (HIRES) The average oxygen 
content of the ^Guzman et al. ( 1996 j ) sources, derived using the Denicolo, Terlevich & Terlevich ( p002) 
expression for 12-|-log(0/H), is 12 -I- log (0/H) = 8.6±0.1. This value is a little bit too high to be completely 



compatible with the majority of the Pilyugin, L. S. (2000) metallicity calculations presented above, although 



it is compatible with some of them. This indicates that the oxygen content is much more likely to lie in 
the upper half of the interval allowed by the oxygen line ratios. It is interesting to note that the HIRES 
measurements of [N II]A6584/Hq; are all around 0.18, and that the 3(T detection threshold at Ha is around 
15% of Ha for the STIS spectra, meaning that the STIS spectra were very close to detect the nitrogen line. 
This is true for both the low resolution and high resolution configurations. 



^In addition, the metallicity gradients of compact systems are likely to be intrinsically small because of the mixing. 



C<1 



in 00 CD 
o o c c 



o o c c 



-H ^ -H -H 

^O'-HiOiNt-C-iCi-i 
0(rocdoOi-lT-H(N 



^ -H S -H -H 
aJt-o~noooicoi>o 

.Sr-id666'-lrHC^ 



00 i> o o 



q q CN rH 
^ 00 00 d d 



o c c 
d d d d 
-H ^ -H -H 

O'^M'-hcOt-hocD 
T-Hxj^rocD'^iCO'-i 

ddddr-lrHc^c^i 



00 l> 



-H 



o o o 

d d d 
^ -1 ^ -H -H 
Ooooi>-^cioo:>oo _ 
S(M_|_|^D':^^oo!r)LO(^^mol^'^l> 
.Sr-lodddwT-HCNcviroc^ooi^ 



^ g ^ 

o . c 

^ d 5 d 

. S S ?8 S 

i d d 9 d 



i d d 



: o _M o 

^ S ^ 3 -H 

o -H o ro 00 00 

S O O . C w 



o ;o ±1 



o 


o 




CO 






o 


o 


O 


d 


d 


d 




-H 


-H 


o 


00 


in 




CO 




d 


d 


d 








q 


q 


q 


d 


d 


d 




-H 


-H 




QO 


CO 






oo 


o 


O 


o 



d --i 



i-H CD CN 

00 00 d 



ID CO l> 00 O O 



lO l> ^ IM 

00 i> d d 



lO C oo l> 00 

(^i (^i 00 i> d 



5 



bo bo bo bo bo bo -.1 ' 

o o o o o o^Z^^c^cq 



3 S 



CI 



.SO 



o a 



j3 g? 5„ 
H O 

CN o 



0) bo 
^> r! n 



n — 

„ 0) .3 

j3 ^ tt 

-0 S B d 



^2 ° 







S 1^ ^ J 



o" 




o" 


_o 




_o 


Oh 




Dh 


ide- 




ide- 








< 






a: 
to 


to 


00 


00 


.lue i; 


value ii 


va 


Si 






[OS 


0) 









■a 2 

6? 



0, ^ 



O ;£3 
CD 



3 



5 ^ 



J3 
ho 





'C 




o 
o 

IN 


o 

T3 

O 


o 
o 

IN 


lund 


idanc 


lund 


Asp] 


abui 


Asp] 




o 

bO 

O 




iinberl 


imberl 






J. 



bo 

B 



J3 



J3 



I 

-d 



O a; 



B 
j3 



bl 



s s 



f3 
T! 



a 

i1 



bo -a 



8 



0) > 



H 5 



^ 2 



s S S 





OJ bl 

.B 



bO 

OJ .cn 



, .B i 

\ V 



g ^ 



K bo o 

B JP f 
« i. 



T3 



ID 



S bl 

(U o 

■05 S 



- 10 - 



3.4. Scaling and Registering of the Data. Calculation of Stellar Colors and NIR 

Luminosities. 

In order to describe the stellar populations found in the galaxies studied, it is needed to match the STIS 
spectroscopic frames with the WFPC-2 images and the NICMOS frames. This is done by forcing that each 
galaxy's centroid in the WFPC2/NICM0S images coincides with the center of the continuum distribution 
of the STIS spectroscopic frames. Once both centroids are calculated, the slit orientation with respect to 
the image is used to determine what pixels in the WF/PC-2 frames correspond to the inner and outer zones 
defined previously for the STIS spectra. 

Figure [2] shows the slit placement in the F814W and F606W WF/PC-2 images. The bidimensional 
spectra presented for objects 1, 3 & 4 are low-resolution G750L frames, whilst the spectrum presented for 
object 2 is a high-resolution G750M frame. In this figure, the spectra have been scaled and approximately 
aligned with the images. It is seen that the STIS spectra only probe the most luminous regions of the 
observed objects, corresponding to the line-emitting region. This area spans less than 1", and it's been 
further divided into two regions. This makes it neccessary the use of the spatial resolution of the HST 
instruments. Even though these objects have been observed from the ground in the UBRIK bands, it 
is not possible to merge those data with the observations presented here because the spatial resolution of 
those observations is not high enough to study the presented objects at the subarcsecond level. The stellar 
population analysis presented in the following section [4] is therefore best suited for the volume within the 
line-emittion regions of the studied galaxies. This region typically has a radius of 1.7 kpc. 

Figure |3] shows the NICMOS fr ames. Each panel is 2.5" on a side. Again, the images have been rotated 
so that the slit silhouette is vertical, as shown in the figure. It is seen that the STIS spectra focus on the 
very central regions of the objects. 

The next step is to determine the B — V color of the areas in the WF-PC-2 F814W and F606W images 
that match the inner (a) and outer (b) zones of the STIS spectra. In order to accurately determine said color, 
several issues are to be taken into account. The first of them is to correct the integrated DN (CCD counts) in 
each area for the imperfections of the charge transfer mechanism. There are several choices to accomplish this 
task. The version used in the WFPC2 Instrument Handbook, based on the 05/31/2002 version of Andrew 
Dolphin's calculations was used. This correction depends on the objects' location in the CCD. It also depends 
on the objects' accumulated DN and CCD gain. This correction is obviously insensitive to the filter used. 
Fortunately, for the observations used, any given galaxy is always imaged in the same CCD location, and the 
accumulated counts are very similar for the two filters used. This correction is therefore very small. It only 
changes the observed colors by less than 0.1 magnitudes, usually much less. The statistical error introduced 
by this correction was estimated by varying the DN parameter of the software tool previously mentioned 
around the observed DN number. It turned out to be that the magnitude error introduced by this correction 
was negligibly small, 0.001 magnitudes in the worst case. For a complete review of this problem and the 



different approaches to correct for it, see Dolphin (2000). 



The second issue is the derivation of the fc-correction. This step is needed to convert from the observed 



I — V colors to rest-frame _B — V^colors. The fits given in Fukugita et al. (1995) for the case of an Irregular 
Magallanic spectrum are used. The fc-corrections are 0.32 for objects 3 & 4, 0.39 for object 2 and 0.47 for 
the furthest object. The errors introduced by this correction are mainly systematic, arising from deviations 
of the real galaxy spectra from the assumed Irr spectrum. These are estimated in around 0.08 magnitudes. 



as was done in Guzman et al ( 1998 1 
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S(kpc) 



3(kpc) 



Fig. 1. — Continuum and emission-line spatial profiles. Only objects 1 & 2 are presented since the spatial 
profiles for objects 3 & 4 were given in Hoyos et al. (2004). The S-axis is placed with respect to the slits 



depicted in Figj2] Increasing S corresponds to upper positions in that figure. The continuum is the thickest, 
solid line; [O III]A5007 is the thick line and the dotted line represents H/3. 



I images. V images. [OH];. 3727. HP.[Oin]U 4959,5007. Ha. 




700A 700A 700A 



Fig. 2. — Matching between the STIS 2-D spectra and the HST/WFPC-2 images. Each spectral strip spans 
2.5"in the spatial direction. The linear scale is given for each object in the figure. Each spectral strip spans 
700A except for object #2, in which the wavelenght coverage is given in each image. The spectral strips 
presented for this objects are not registered with respect to the spectral strips shown for the other objects 
for this reason. The Ha strip for object 1 is not perfectly registered with respect to the Ha strips for objects 
3 & 4 because the Ha line fell very near to the CCD edge. This spectral strip still spans 700A, though. The 
slit orientation is shown superimposed on HST F606W and F814W images which are 2.5"x2.5"on a side. 
It is interesting to note that the spatial distribution of the continuum is lopsided with respect to the peak. 
This is most evident in objects 3 & 4. 



-VI- 




Fig. 3. — NICMOS F160W reduced images. Each postage stamp is 2.5" on a side. The sht is placed as in 
Figure |2] The number in the upper left corner of each panel corresponds to the identifying numeral of Table 

m 
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The third issue is to derive the emission-line free colors. The method presented in Telles & Terlevich 



(1997) is adopted. This method uses the rest-frame equivalent widths, together with the B and V filter 
transfer functions to estimate the contribution of the emission lines to the observed galaxy colors. The 
method is now sketched. Given any observed color A ~ B, it can be written as: 



B-A = Cba- 2.5 xlog^^±^ (1) 

+ ^A 

where Cba is the photometric constant, L'^ is the continuum luminosity in the X band and is the line 
luminosity falling within the X band. 

If the line luminosities are considered as small perturbations to the spectrum, the above equation can 
be cast into the following form. 

B-A = {B-A)o + 1.086 x {L'JL'^^ - L's/L%) (2) 
where (B — A)q is the color that would be observed in the absence of the emission lines. 
It is also possible to write the color correction as: 

A- B = (A- B)o + l.086x ( — ^ ^ — ^ (3) 

^ ' \IaxWa IbxWbJ ^' 

where Wx is the bandpass of the X band, and Ix is the average specific intensity. 
The color correction can therefore be written as: 



1.086 X (L^^ _ EWl) \ 

\ Wa Wb J ^ ' 

where EW(i) is the line equivalent width of the i*^ emission line, and Xi is the filter transfer function 
evaluated at that wavelength. 

It is important to note that, even if the line equivalent widths used in the above formulae bear large 
uncertainties, any error in the measured equivalent widths is greatly diminished or cancelled out. This 
happens because the same equivalent width appears in both terms of the color correction with opposite 
signs, although with possibly different multiplicative factors. 

The applied corrections range from essentially zero to 0.25 magnitudes. The statistical error introduced 
by this correction is very small. Given the derived uncertainties in the measured equivalent widths, these 
errors are less than 0.02 magnitudes. 

Finally, the ^{B — V) color excess is derived from the c(H/3) values for each zone. These color excesses 
are very small and are virtually error free. The resulting colors are therefore rest-frame, stellar population 
only, extinction corrected i3 — V^colors. The total photometric errors, obtained adding in quadrature all the 
aforementioned uncertainties, are estimated in 0.09, for all objects. These colors are gathered in Table [2] for 
the inner and outer parts of the line-emitting region in which the stellar populations are studied. 

The NICMOS data were used to calculate the total F160W luminosities of the line emitting regions. It 
is needed to use the NICMOS data in this way because of the larger angular extent of the NICMOS Point 
Spread Function, that makes it impossible to calculate the observed fluxes of the (a) and (b) zones separately 
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with the same accuracy as for the optical frames. The measured fluxes and luminosities are gathered in Table 
|2] The reported luminosities all assume c{H[3) = 0.2, and they also assume that Ax/ Ay = 0.4, where Ax 
is the extinction in the wavelength range that maps into the F160W passband. This value is the average 
value for that ratio between the / and J bands. The contribution from the nebular continuum, which begins 
to be important at these wavelengths has been estimated at less than or around 10% using the standard 



techniques presented in Osterbrock (19891 



It is seen in Table |2] that the colors of the inner zones are usually bluer than those of the outer zones of 



the line emitting region, in accordance with the Guzman et al ( |1998^ color profiles. This can be interpreted 
in the light of the results presented in subsection |3.2[ Whenever one region in a galaxy is redder than other 
region, it can be caused by differences in reddening, metallicity, or age. The first of these factors is unlikely 
to play a role in this particular case since the reddening has been estimated, and it is very small. The second 
factor is also unlikely to be responsible of the observed color gradients since the metallicity differences are 
likely to be very small in these systems at such short spatial scales because of very efficient mixing by 
supernova driven winds. It is moreover concluded that the observed color differences between the (a) and 
(b) zones are more likely to arise from age differences between these two zones. Specifically, the underlying 
stellar populations in the (b) zones are typically older than the stellar populations of the (a) zones of their 
respective galaxies. The chemical enrichment histories of both regions have probably been very similar for 
the reasons outlined above. This is an important finding for these distant and compact sources. 

Alternatively, it might be possible to think that the color gradient is caused by mass segregation. If , 
for instance, most of the newly formed ionizing stars were in the (a) zone, the color of the inner part would 
be naturally bluer than the color of the outer zone. This scenario could also explain the observed differences 
in the ionization ratio between the two zones, which would be higher in the inner zone because of the higher 
amount of ionizing photons available to ionize the inner parts of the nebula. Although mass segregation is a 



concern for massive clusters (see, for example Stolte et al. (20051 and de Grijs et al. (2002)), it only operates 



at very short distances (10 pc at most), and not at the much larger scales of interest in the present case. 
This possibility of a spatially varying Present-Day Mass Function can then be safely ruled out. 



4. The Evolutionary Models. 



The observational data gathered for the LCBGs presented here can be used to construct galaxy models 
that represent the objects. Although the models presented are very crude and tailored, they provide at 
least an idea of the stellar inventory and possible evolution of these LCBGs. This will shed light on the 
possible evolution and age of these galaxies. In this section, the star formation history of the four LCBGs is 
studied using evolutionary population synthesis techniques. The models have been built using the ,BruzuaIl 
& Chariot (20031 and Leitherer et al. (19991 evolutionary models. For each galaxy and extraction, (a) 
or (b), two simple stellar populations are combined. One population is designed to mimick the observed 
properties of the ionizing stellar population, and the other population is chosen to match the observed colors 
when combined with the ionizing population. This latter population therefore represents the underlying 
population. We note that the population analysis which will be used here is carried out for both the inner 
and outer zones of the line-emitting region separately, taking advantage of the HST spatial resolution. In 
addition, the fact that the star forming region is split into two subregions can also help to model possible 
delays between the star formation episodes in the inner zone and in the outer zone of the line-emitting region. 



The metal content of the ionizing stellar population is assumed to be equal to the metal content of the 
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ionized gas. Since the only information available on the metal content of the ionizing population of these 
sources comes from the gas phase oxygen abundance^ two sets of models are constructed. Models #1 
assume that the ionizing population has a metallicity oi Z — 0.008. Models #2 adopt Z — 0.004 for the 
metal content of the ionizing population. These two numbers can be mapped into 12 + log (0/H) = 8.4 and 
12 + log(0/H) = 8.0, approximately, assuming solar chemical proportions. Therefore, for each galaxy and 
extraction, two models of different metallicities are built. The stellar inventory of the observed galaxies is 
described according the following steps: 



1. The first step consists on the derivation of the masses of the ionizing populations. 

2. The second step is the determination of the ages and colors of the ionizing populations. In this step 



the Leitherer et al. (19991 models are used 



3. In the third step, the mass of the underlying stellar population is derived, together with its age. The 



Bruzual & Chariot ( 2003 1 models are used to represent the older stellar generation. 



4. The last step is a double check of the derived masses using the measured NICMOS luminosities. 



4.1. The mass of the ionizing populations. 



The mass of the ionizing population is a very important parameter that defines it. If the mass of the 
ionizing generation is high enough, the initial mass function will be well sampled. If the initial mass function 
is not well sampled, it might be the case that the models can not be used to describe the stellar population. 
This effect is commonly known as the richness effect or also as the poisson shot noise. Depending on the 
slope and cutoff masses of the IMF used, this can become a problem for population masses lower than 
1O^M0, but for the usual Salpeter IMF the mass limit is lower. This issue has been extensively studied in 



the series of papers Cervino et al. (2000 2001 20021 



The mass of the ionizing population, M.^ can, in principle, be derived from the luminosities of the 
observed recombination lines according to the following considerations. As the stellar population ages, the 
number of hydrogen ionizing photons per unit mass of the ionizing population also decreases. Assuming 
that the age of the ionizing population is related to the intrinsic W^p, a relation should exist between the 
number of hydrogen ionizing photons per unit mass of the single stellar generation and Wp. Such relation 



is given for single- burst models in Di'az et al. ( 2000a I . The expression used is 



log{Q{H)/M^) = 44.8 + 0.86 x log W^, 



HP 



(5) 



where Q(H) is the number of hydrogen ionizing photons per second, M* is expressed in solar masses, and 
Whp is the "intrinsic" equivalent width, the one that would be observed if no underlying stellar population 
was present. Again, this expression works best in stellar populations with well sampled IMFs. The number 
of ionizing photons is an integral over the mass range more affected by the richness effect, and it is therefore 
very sensitive to small number statistics when applied to low mass populations. The IMF used to derive this 
expression had an exponent of 2.35, a a lower mass limit of O.85M0 and a upper mass limit of 120 M0. 



'^12+log(0/H) lies between 7.8 and 8.4, although the preferred value is 8.4 because of the 
measurements. 



Guzman et al. 



(19961 HIRES 



- 16 - 



However, the presence of the underlying stellar population introduces a degeneracy in the measurements 
of the equivalent widths. This prevents us from measuring directly the equivalent widths that would be 
observed if only the ionizing stellar population was present. The main goal of the stellar population analysis 
presented here is then to unravel the different stellar populations. This is carried out following the ideas 



presented in Di'az et al. (2000b). In order to solve this degeneracy in equivalent width introduced by 
the presence of the underlying population, the upper envelope of the relationship between log VFh/3 ^^nd 



log([OIII]/[OII]) presented in ( jHoyos fc Dfaz|[2006| Fi gure 6 of) was used to estimate the equivalent width 



that would be observed in the absence of any previous stellar generation. This "effective" equivalent width 
is called W^. This assumption is now briefly justified. In the Wh/3-[0 III]AA4959 + 5007/[O II]A3727 ratio 
diagram for very luminous H II galaxies, if the observed range of ionization degrees is ascribed to the age of 
the ionizing population, it is reasonable to assume that the vertical scatter shown by the data is caused by 
different contributions of continuum light from the underlying populatioii|^ The expression used is: 

log W^fj = 2.00 ± 0.02 + (0.703 ± 0.028) x log [OIII]/[OII] (6) 

The upper envelope of this relation, which defines the "effective" equivalent width, is then relevant to 
estimate the "intrinsic" equivalent width because the objects located on top of this upper envelope are likely 
to be the ones in which the underlying population contribution is minimum. This "effective" equivalent 
width allows us to estimate the mass of the ionizing population assuming that is equal to the real or 
"intrinsic" equivalent width that would be measured if the starburst could be observed in isolation. This is 
done via the expression given above using the measured H/3 luminosities to obtain the number of hydrogen 



ionizing photons per second through standard nebular analysis techniques Osterbrock (19891. The resulting 
stellar masses of the ionizing population are given in Table [3] Typical ionizing population masses range from 
8 to 30 million solar masses. The systematic errors affecting these calculations are explained in the following 
paragraphs. 

The "effective" equivalent widths, which are mostly in the lOOA- 200Ainterval are still lower limits only 
to the real value of the equivalent width because the empirical method used here to derive the former does 
not forbid the existence of objects with higher observed equivalent widths for their observed ionization ratio 



[O III]/[0 II]. This is indeed a possibility since the work presented in Hoyos & Di'az (2006 1 used local galaxies 
and the galaxies analyzed here are intermediate- reds/ii/f sources. Therefore, the derived value for the mass of 
the ionizing populations obtained using Q{H) (or, rather, the H/3 luminosity) and can only be only lower 
limits to the real masses. It has to be said, however, that the possible differences between the "intrinsic" and 
the derived "effective" equivalent widths can not be very high. The maximum attainable equivalent width, 
which corresponds to the case in which only the nebular continuum is considered, is lOOOA. When the stellar 
continuum is taken into account, the maximum equivalent width then drops to about 400A, depending on 
the IMF, metallicity and age of the population, and the largest H(3 equivalent widths ever measured are 
close to this latter quantity. This implies that the "intrinsic" equivalent widths of the observed starbursts 
have to be in the 200A- 400Arange. These considerations allow us to estimate the possible systematic errors 
introduced by this approximation. Using the typical "effective" equivalent width, and an average value for 



^There is, however, a metallicity effect. In the absence of any underlying stellar population, H II galaxies with higher 
metalicity clouds will present a lower ionization ratio due to lower effective temperatures of their ionizing stars. However, there 
will be no change in the equivalent width of H/3, to zero-th order. Therefore, should be a function of both the ionization 
ratio and metallicity. Only the ionization ratio dependence of the "effective" equivalent width will be used here. 
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the maximum "intrinsic" equivalent width, the real masses of the ionizing populations could be up to a 
factor of 2 larger. However, even if the "intrinsic" value could be observed or somehow derived, the mass 
of the ionizing component obtained would always remain a lower limit because some photons could actually 
be escaping from the nebula, and dust globules will always be present amidst the gas clouds. Moreover, if 
the ionizing population turned out to emit more photons harder than 1 Rydberg than the number predicted 
by the models used to derive Eq. |5j the mass would be underestimated, too. The new Geneva models with 
rotation at low metallicity ( Maeder fc Meynet|2001 1 point towards this direction. For these reasons, further 
steps towards the determination of the "intrinsic" equivalent width from the "effective" one were not taken 
and both quantities are assumed to be equal. 



4.2. Determination of the Age and Intrinsic Color of the Ionizing Population. 

The value of the "effective" equivalent width can be used to obtain more information about the under- 



lying stellar populations via the Leitherer et al. (1999) models, making some assumptions on the metallicity 



and IMF of this stellar generation. Again, the "effective" equivalent width allows us to estimate, or at least 
to place significant lower limits to, the age and intrinsic color of the ionizing population assuming that Wj^ is 
equal to the real or "intrinsic" equivalent width that would be measured if the starburst could be observed 



in isolation. The Leitherer et al. ( 1999 1 models are used now 



Figure |4] presents nine models from the Starburst99 ( Leitherer et al.||1999 l library, which represent the 



evolution in the log vs.B — Vplane of a single stellar population from its birth till it is lOMyr old. The 
nebular continuum emission is included in these diagrams. This evolution is the solid line in those plots. 
Ages increase from top to bottom. Three different metallicities and IMFs are presented, as indicated in 
the figure. The derived for both the inner and outer zones of the line-emitting regions of the observed 
sources is plotted against the emission-line free, reddening- corrected color in the case of the studied LCBGs. 
This is the color that would be measured if the stars could be observed in isolation. 

It is interesting to note that the values of the "effective" equivalent width are close to the maximum 
equivalent width value of the plotted curves, which is the maximum value of the "intrinsic" equivalent width 
for zero age clusters of the corresponding metallicity and IMl|^ further highlighting the fact that, even 
though the "effective" equivalent width is a lower limit to the real value, it is a reasonable approximation to 
it, specially if one considers that, even if the real value could be somehow known, it would still yield lower 
limits to the mass of the ionizing clusters because of photon escape and dust issues. It is clearly seen that 
3/4 of the regions studied require the presence of an underlying stellar population since they do not lie along 
any solid line. They can not be reproduced by a single stellar population, then. It is also seen that the 
single stellar populations with a truncated Salpeter IMF have maximum "intrinsic" equivalent widths which 
are lower than the derived values of Wp of many regions. This means that we can rule out these models as 
descriptions of the ionizing stellar populations found in the observed LCBGs since the "effective" equivalent 
widths are in any case only a lower limit to the real one. This can be interpreted in the sense that the 
observed galaxies have ionizing populations so massive that they fully sample their respective IMFs, which 
is a necessary condition to achieve high equivalent widths. In addition, the models whose IMF has a slope 
of 3.33 are redder than the inner region of object 3, which makes it impossible to accomodate any older and 
redder underlying population in this model and there are also some regions with greater than the zero 



About 350A, for the presented cases. 
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age equivalent width of those models. The models with Salpeter IMFs are therefore chosen to describe the 
ionizing populations, simply because they are compatible with all these restrictions. Although the reasons 
to discard the m^^-^^ IMF are of course weaker, considering multiple IMFs is by no means justified given 
the few inputs data at our disposal. The use of the Salpeter IMF provides a crude yet consistent description 
of the ionizing population found in the observed regions. 

Assuming that Wp equals the real equivalent width, as was done before, and using the Salpeter IMF 
for the two metallicities that correspond to the mctallicitics of the two models #1 and #2 used, it is 
possible to calculate the ages and intrinsic colors of the ionizing population. They are given in Table |3] 
The ionizing clusters are found to be very young (less than 5Myr in most cases). As before, the use of 
the "effective" equivalent width instead of the unknown "intrinsic" value introduces a systematic error in 
the age determinations. The ionizing populations are then bound to be even younger than the presented 
quantities. When dealing with very young coeval single stellar populations, it has to be kept in mind that 
the observed properties of such populations do not change to a great extent until they are about 3Myr old, 
simply because of the fact that this is the minimum lifetime of even the most massive stars and also because 
the approximation of an instantaneous burst might fail at such short timescales. The main conclusion is 
that the ionizing populations of the observed regions are probably extremely young. The systematic error 
in the intrinsic color of the ionizing population is also estimated from the color scatter of the models used 
between zero age and the maximum age allowed by the observed Wa. All errors are 0.04 magnitudes. 



4.3. The Age and Mass of the Underlying Population. 



The next step to describe the stellar populations found in the observed LCBGs is then to calculate the 
age and initial mass of the underlying stellar population. The single stellar populations used to describe the 



underlying generation of stars are chosen from the Bruzual & Chariot ( 2003 1 library. The evolutionary tracks 



used are the Padova 1994 ones and the IMF chosen is the Salpeter IMF. The metallicity of this underlying 
stellar population is Z — 0.004 for the case of models 1 and Z — 0.0004 for the case of models 2. These 
are lower than the metallicities used to describe the ionizing populations. In this step, the ionizing cluster 



is represented by the single stellar population available in the Bruzual & Chariot (20031 library that most 

( 1999 ) model used before in terms of metallicity, IMF and W^-derived 



Leitherer et al. 



closely resembles the 
age. 

Assuming that the metallicities of the ionizing populations are higher than the metallicities of the 
underlying populations tries to emulate the chemical evolution of the observed galaxies since their inception 
till the onset of the observed starbursts. Although the choice adopted is not justified by the use of any 



particular chemical evolution model, it is at least a zero-th order solution to the problem since the Bruzual 



& Chariot (2003 1 models only offer a limited array of metallicities for the different single stellar populations. 



In the framework of the closed box model and the instantaneous recycling approximation, the corresponding 
gas consumption fractions for the required metal enrichments are 10% for models #1 and 60% for models #2, 
assuming zero metallicity in the distant past, before the starbursts that gave birth to the underlying stellar 
generations took place. The gas consumption fraction is larger in the metal poor models since the difference 



in metallicity between the two simple populations mixed to construct them {Zio^ = 0.004, Z^ 



idcil 



0.0004) 



is much larger than the metallicity difference of the two simple populations combined for the first models 
(Zion = 0.008, Zundori = 0.004). The yields provided by [Molla | ([2006[) and [Maeder] (|1992|) were used in 



these latter calculations. Although these estimates, by themselves, do not prove that this picture is accurate 
to describe the star formation history of these objects, they show that the assumed metallicity differences 
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Fig. 4. — Evolution of a single stellar population in the logWHi3-B — V plane from t = OMyr to t = lOMyr. 



The data were taken from the Leitherer et al. ( 1999j ) library. In all cases, an instantaneous burst is assumed, 
and the nebular continuum emission is taken into account. Two different oxygen contents {Z — 0.008, models 
#1, lower panels, and Z — 0.004, models #2, upper panels) and three different IMFs (Salpeter from I.OM0 
to IOOM0, left panels, a power law with an exponent a — 3.3 with the same mass limits, central panels, and 
a Salpeter IMF truncated at 30 Mq, right panels) are plotted. For the LCBG points, is plotted against 
the emission-line free, reddening corrected color. The horizontal error bars represent the 0.09 magnitude 
uncertainty, and the vertical semi error bar reaches till the observed equivalent width, in order to indicate 
the importance of the contribution of the underlying stellar population to the optical light. The big black 
dots represent the 0,2,4,6,8,10x10^ years checkpoints. Object 1 is represented by half-filled squares, object 
2 is represented by 0> the location of object 3 is given by pentagons and x denotes object 4. Two points 
are given for each object, corresponding to the (a) and (b) extractions. In all cases, the (b) apertures have 
redder colors. 
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between the stellar populations used to build thew models are, at least, possible. In addition, it has to be 
kept in mind that the global remaining gas fraction is bound to be lower for the metal rich models in the 
framework of the closed box model and the instantaneous recycling approximation than for the metal poor 
ones. In particular, given the assumed metallicity for the ionizing population of the second set of models, 
the gas consumption fractions of pristine gass have had to be higher than 60% in both cases, again in the 
closed box with instantaneous recycling approximation. 

The age and mass of the second stellar population are then chosen so that the predicted equivalent 
width by the time the current starburst is observed matches the measured equivalent width and the total 
color equals the calculated rest-frame, stellar population only, extinction corrected i? — V^color. In essence, 
this procedure estimates the amount of "red light" that has to be added to the "blue light" of the ionizing 
population in order to match the two observables available (iJ — ycolor and measured H/3 equivalent width) 
at the time in which the observed starburst takes place. This second stellar population is therefore designed 
to match the older stellar population. 

The difference between the observed B — VcoIot and the previously derived intrinsic i? — X^color of the 
ionizing population (given in Table [s]) can be approximately matched with "how red" , and the difference 
between the observed equivalent width and the "intrinsic" equivalent width roughly translates into "how 
much" red light has to be added, although both quantities are of course interlaced. 

The errors in the observed equivalent widths (typically 17%) and stellar colors (0.09 magnitudes, as 
estimated above). Translate themselves into errors in the relative mass of the underlying stellar population 
with respect to the ionizing population and age. Although it is not true 100% of the error in the relative 
mass comes from the uncertainties in the measured equivalent width and that 100% of the errors in the age 
of the underlying stellar population arise from the errors in the derived stellar colorj^ these are the main 
dependencies. Assuming this simplified error scheme to hold, the typical uncertainties in the relative mass 
of the ionizing cluster with respect to the mass of the ionizing cluster are 20%. The uncertainties in the age 
of the underlying population are 25% for ages below lO^yr, and twice this amount for older ages. This is so 



because of the flattening in the B — V-a,ge relationship of the Bruzual & Chariot (20031 models. The results 
are also presented in table |3] 



4.4. NICMOS NIR Photometric Masses. 

The final step in this crude description of the stellar populations found in the observed LCBGs is made 
possible by the use of the NIR data. This information is now employed to check the likelihood of the derived 
stellar populations, by using the total luminosities emerging from the line-emitting regions in the wavelength 
range that transforms into the F160W wavelength interval measured by NICMOS to rederive the total stellar 
masses contained in the line-emitting region, which is the area probed by the STIS spectra. This is done 
by calculating the relevant M/Lx ratios via the M/lJ^ ratios, making use of the programs presented in 



Bruzual & Chariot (2003). The statistical uncertainty that affects these NIR derived photometric masses is 



a combination of two effects. The first effect is the intrinsic uncertainty in the flux measurement, and the 



^"In the sense that, for instance, iso-agc contours in the age-relative mass plane do not map directly into iso-W^ contours in 
the W f}- stellar color plane. 

^^The Lg/Lx ratios are calculated for the (previously derived) model that best describes each particular object. The two 
families of models #1 and #2 are used. 
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second effect comes from the scatter in the L b / Lx ratio. This latter error source dominates over the first 
one, having a scatter of 30%. This number was derived by recalculating the Lb/Lx number for different 
BC03 spectra of the two metallicities considered for the underlying stellar populations, with the typical ages 
predicted in the models previously constructed. 

Total model masses, NIR photometric masses as well as the ratio between the two are given in Table |3] 
It is seen in that table that the first set of models produces NIR photometric masses that are in very well 
agreement (to within 25%) with the masses derived using only the spectroscopic and optical photometry. 
The agreement for the second set of models is somewhat worse (to within a factor of 2) . These results lead 
to think that the models ^1 are a better overall description. 



4.5. Mental model for the observed objects. 

Table |3] lists the ages and masses of the stellar populations that best represent the observed objects in 
the framework of the presented models. This table also gathers the star formation rates and the percentage 
of the dynamical mass that is accounted for by the stellar populations used. 



The star formation rate (SFR) was calculated from the Ha luminosity as in Kennicutt (1994), valid 
for Tg = 10"^ K and case B recombination (all the ionizing photons are processed by the nebular gas). The 
results are given in Table [3j 

SFR{MQyr~^) = 7.9 x 10"'*^Ljj„(er<? s"^) (7) 

Although the preceeding SFR estimate was based on spatially averaged SFR for disk galaxies, we believe 
it can be applied for the current sample of starburst galaxies. The Ha derived SFR only probe the most 
recent star formation episodes of any given galaxy because the Ha line becomes very faint soon after the 
more massive stars, born in the more massive starbursts, die out, and this is true for all normal galaxies. 



For this reason, it is possible to use the Kennicutt (19941 expression for the global SFR of actively star 



forming galaxies, too. In fact, the SFR Ha estimate has been calibrated against many other instantaneous 



star formation tracers (Rosa-Gonzalez et al. [2002 ). Differences between the Kennicutt ( 1994( ) Ha-derived 



SFR and the real SFR of the studied galaxies are bound to arise mostly from other issues, such as differences 
in the IMFs used. 

The calculated star formation rates are in the range from 0.8 to 6.8 MQyr~^, indicating that the star- 
forming episode is very strong (30-Dor SFR is 0.1 Moyr~^). It has, however, to be said that the line-emitting 
regions typically span around 3 kpc, and the 30-Dor line-emitting region is a only few hundred parsecs. The 
star formation rates per unit area are therefore of the same order of magnitude. 

Tables [3] shows that the ionizing populations are very young, less than 5Myr in almost all cases, although 
the average age is slightly over 3Myr. The ionizing clusters are then largely unevolved. The resulting ages in 
both model sets indicate that these galaxies are being observed at the peak (or a little bit past of) of their 
line luminosities. It is also the case that the ages of ionizing populations are very similar for the inner and 
outer regions of each galaxy. 

The inferred underlying populations are very different amongst the studied galaxies. There are systems 
with fairly young non-ionizing populations (around 300Myr old), but there are LCBGs whose underlying 
populations are very old (older than l.OGyr). In the case of the inner region of SA57-10601, the estimated 
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stellar population is compatible with a single ionizing stellar population. The modelled age of the older 
stellar population is usually larger for the outer regions than for the inner regions of each object. However, 
given the uncertainties in the ages estimates for the underlying populations (between 25% and 50%), it is 
not possible to extract much information from these ages alone. The average age of the underlying stellar 
populations found are 700Myr for both models #1 and #2. 

The derived stellar masses account for 20% to 60% (the average value being 40%) of the virial masses 
within Re, according to the first set of models. In the case of models 2, the stellar mass accounts for 4%- 
70% (the average value is again 40%) of the derived virial masses, where the virial masses were determined 
as in Hoyos et al. (20041, assuming the measured gas-phase velocity dispersion equals the star's velocity 
dispersion and that LCBGs are in dynamical equilibrium. These are the virial masses within Re- Since 
the mass estimates yielded by the models used here are only lower limits for the reasons outlined abovep^ 
there is not much room left for the presence of dark matter in the innermost regions of these objects. The 
dark matter would, in any case, be spread out in a vast halo around the objects by its very nature. The 
presence or not of significant amounts of dark matter in the optical cores of these systems is important. If 
the conjecture introduced here that these systems are not dominated by dark matter in their very centers 
turns out to be true, it would imply that these systems could not possibly evolve into today's population 



of Low Surface Brightness objects. A very simple and crude estimate based in the Mamon & Lokas (20051 



models and the Navarro et al. (20041 simulations indicate that the mass of dark matter within a 1.5 kpc 
spherij^ is 2 x 10^ ±40%, assuming that the total mass of the dark matter halo is ten times the measured 
kinematic mass. This contribution is much less than the kinematical masses derived from velocity dispersions. 
In order to determine the dark matter content and density profile of these sources with greater accuracy, 
measurements of the line of sight velocity distribution out to several (around 5-9) effective radii would be 
required. Such data are of course not available for the distant objects studied here. 

It might be tempting to think that the remaining 60% of the dynamical mass not accounted for by the 
stellar component might be hinting the existence of a very old stellar population, too faint and red to leave 
a trace in the B — ycolor used here but as massive as the already detected underlying stellar population. 
However, this hypothetical phantom population can not be older than ~ 5Gyr, simply because of look-back 
time reasons, implying that such very old and massive population can not possibly exist. Such population 
would still be shining in the _B — l/coloi[^ because the turn-off point of this population would be around 
solar- type stars or slightly earlier. More strongly, the good agreement between the optically derived masses 
and the NIR derived masses sets very strict limits to the luminosity of this possible population, virtually 
forbidding it for the case of models #1 and for some of the models #2 as well. Finally, the presence of 
significant amounts of neutral or molecular gas in the central regions of the sources studied can be safely 
discarded. It is very likely that the extreme starburst taking place at the cores of these objects is well capable 
of fully ionizing the ISM in the innermost zones of these galaxies. 

It is therefore concluded that the stellar inventory found using the simple models devised here can not 
be very far from reality and that the disagreement between the stellar masses and the dynamical masses 



^ We remind the reader that the masses of the underlying populations are scaled with respect to the mass of the ionizing 
population. This is calculated via W^, which is in turn a lower limit to the real equivalent width. Other factors such as photon 
escape, dust and enhanced photon production have an impact, too. 

^■^ Similar in size to R^. 

Assuming these stars occupy a volume similar to the other populations, i.e., the surface brightness of this population is not 
lowered for other reasons 
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arises from the fact that the mass calculations for the stellar component are forcibly lower limits. 



The models constructed here are similar to the two-bursts models presented in Guzman et at (19981. 



In the models presented here, the ionizing cluster is less massive than the ionizing clusters used in that 



work. In addition, the underlying stellar populations used in Guzman et at ( 1998 ) are much older than the 



ones derived using the models presented here. The relative masses of the ionizing populations considered in 



Guzman et al ( 1998 1 went from 1% to 10%, while the results presented in this paper lead towards much lower 



relative masses for the newer stellar generation, between 0.1% and 3%. Another very important difference is 



that the underlying stellar population used in Guzman et al ( 1998 1 was set to be 3Gyr old at the time of the 



ionizing generation birth, whereas in the current work this age is found to be much shorter, around IGyr. 
Finally, this section is closed with a brief summary of how the models are built. This procedure is 



carried out for each region (a) and (b), and assumed metallicity of each object. See Diaz et al. (2000b I for a 
complete account on how to mix several single stellar populations to match the observed properties of star 
forming regions. 



• The first step is to derive the II/3 "effective" equivalent width. This is done via the ionization ratio 
[O III]/[0 II] and the relationship given above. 

• This "effective" equivalent is assumed to be equal to the "intrinsic" value and used to determine the 
age and i? — ycolor of the ionizing stellar population. In this step the error budget is dominated by 
the systematic error of assuming that both quantities are equal. 

• Using the measured H(3 luminosity, the mass of the ionizing population is estimated. This estimate 
can only be a lower limit to the real value because some unprocessed photons might be escaping from 
the nebula, and dust globules will be present. This also translates into the derived masses for the 
underlying stellar populations. 

• Using the emission-line free _B — V^color and the observed Hj3 equivalent width, the age and relative 
mass of the noionizing population are derived. The errors in these quantities are dominated by the 
uncertainties in the required color and required specific intensity at 4860Aof this population. 

• The NIR observations are then used to check if the models created to describe the optical properties of 
the observed galaxies also describe their NIR properties. It is the case that the agreement is excellent 
for the metal-rich models, and somewhat worse for the metal-poor ones. 



5. Comparison with Local Galaxies: Bright, Local HII Galaxies as Nearby LCBGs and the 

LCBG-Spheroidal Connection. 

It is interesting to compare the evolution of the models described here with present-day galaxy samples. 
In this section, several predicted observables are compared to the properties found in local systems. In 
particular, the issue of whether or not the stellar populations predicted by the presented models are similar 
to the underlying stellar populations of local HII galaxies or in the less luminous Blue Compact Dwarf 
Galaxies is addressed. This section also deals with the passive evolution of LCBGs, and their connection 
with the local population of dwarf elliptical systems. 



The HII galaxy sample used here was taken from the Telles & Terlevich (19971 sample. This sample 
comprises 15 galaxies, which were selected from the Spectrophotometric Catalogue of HII galaxies (Terlevich 



et at. |l991 SCHG). Their equivalent widths are between 26A and 170 A and their Ms ranges between -14.5 
and -21.5, although with a preference for the brighter values. The measured velocity dispersions range from 
16 km s^^ to 50 km s~^. The galaxies also exhibit a wide morphological variety. This sample was chosen 
since HI I galaxies also show conspicuous and galaxy- wide emission line spectra. They are quite similar to 



the LCBGs presented here in this respect, as hinted by Koo et at. (19941, Guzman et al. (19971 and Hoyos 



|et al.| ( 2004| . Their small radii and velocity dispersions indicate that they ought to be as massive as distant 
LCBGs. In addition, they are also known to posses violent star forming systems, just like intermediate- 2; 
LCBGs. It is then natural to wonder if bright, local HII galaxies are nearby copies of distant LCBGs or if 
they are different from the distant sources in some respect. 



The BCDG control sample is the one found in the Doublier et al. (1997 1999) works. This sample 



consists on 44 galaxies, which were selected from the Second Byurakan Survey ( Markarian [1967 1 and other 
lists. The Caon et al. (2005) BCDG sample was also included. It adds 7 objects. Their luminosity profiles 



are typically r^^"* profiles, although there is a sizeable fraction of purely exponential profiles and mixtures of 
both profiles. Their star forming regions are clearly detected in B — R color maps and show an age spread. 
The presence of very old red giant stars indicates that these systems began to form their stars several Gyr 
ago, at the time in which the intermediate redshift LCBGs were observed. This sample was chosen since 
BCDGs could be considered to be the low-luminosity end of HII galaxies in many respects. Therefore, they 
can be used to make educated guesses on the properties of other populations of faint intermediate redshift 
objects whose luminosities are not high enough to be considered LCBGs and are therefore very difficult to 
study observationally, given their distances. 

It is also possible to investigate the passive evolution of these systems in the framework of the presented 
models. This scenario of a passive evolution assumes that no further episodes of star formation takes place 
from their look-back times to the present day. Given the high intensity of the observed starburst, it is 
not unlikely that the interstellar medium will be cast out of the galaxy, preventing further starbursts from 
happening. However, it has to be borne in mind that the look-back times of the LCBG sample used here are 
^ 5 X lO^yr. This means there is a huge amount of time for they to undergo other star formation episodes, 
mergers or other phenomena that would, of course, not be represented by the simple passive evolution 
scenario used. In particular, the evolution of the intermediate redshift LCBGs studied here into today's 



dwarf elliptical systems is studied. The dE/Sph sample used for comparison purposes is that of jCuzman 



(1994). We also use four other well-known dwarf elliptical galaxies from the local Universe. The observed 



LCBGs galaxies are also compared with The Large Magallanic Cloud. 

It is important to note that the available colors and magnitudes for the comparison samples are integrated 
data while the models for the LCBG sample have been constructed from the line ratios and colors of the 
line emitting region of each object. The model results for the inner and outer regions of each galaxy have 
been therefore mixed together and scaled using the appropriate weights to derive the properties of the 
whole sources. The main assumption is that the underlying stellar population of the central kpc of these 
objects is qualitatively similar to the stellar population found in the outskirttj^of each galaxy. The absolute 
luminosities were further corrected by multiplying them by a factor of two, to account for the luminosity of 
the stars found in the outskirts of the studied object This way, the observed galaxies can be represented 
by a single point each in color-color or color-magnitude diagrams. 



Here, the word "outskirts" means "outside the Hne-emitting region." 
This factor is needed because the STIS spectra only samples within Re 
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Figure |5] presents the evolution of the stellar populations described by the models used here from their 
predicted creation time to the present day in the Mb-B — R plane for the two sets of models used. The 
B — R color information is an output from the .Bruzual & Chariot (20031 programs. The B magnitude is 



the standard Johnson magnitude while the R magnitude is the Cousin's R magnitude. Each intermediate 
redshift LCBG is represented by a colored line resulting from appropriately combining and scaling the models 
for the inner and outer zones of the line-emitting region, as stated above. This is also true for the rest of the 
diagrams. The Sph sample is represented by simple filled dots and the local H II galaxies are represented by 
half-filled circles. The BCDGs from the Doublier et al. (1997 19991 samples are depicted as crosses. The 



empty squares represent the Caon et al] ( [2005 ) galaxies. 

The B — R colors plotted in figure [5] for the H II galaxies and BCDGs are not the B — R colors of the 
whole galaxies, rather, they are the B — R color of the underlying stellar component, excluding the starburst. 
This B — R color has been estimated as the terminal color in B ~ R color profiles, or from the V — R colors 
of the extensions, assuming _B — 1^=0.65. The B — R color of the dE sample is the color of the entire galaxy, 
since there is no ionizing population in these objects. Therefore, it is possible to compare both B — R color 
estimates in the sense that, if the B — R color of the wings of an H II galaxy or BCDG is similar to the 
color of a Sph system, it would be possible to think that the starbursting H II galaxy or BCDG is a Sph 
system experiencing one further starburst. Figure |5] shows that the underlying stellar population of H II 
galaxies and most BCDGs are younger than the average B — R color of Sph systems, implying that there 
are have to be other differences between these galaxy types. In the case of the H II galaxies and BCDGs, 
the plotted values of Mb are the blue absolute magnitudes of the stellar component only. These values were 
derived from the global Mb using the appropriate available data in the references used, and were corrected 
for galactic extinction according to'Burstein & Heiles (19821. The values of Mb for the LCBG models were 



normalized to the starburst's luminosities, which are common to both the metal rich and metal poor models. 



The distance modulus used to derive the blue absolute magnitudes of the Guzman (19941 sample from 



their reported apparent magnitudes is 35.0 ± 0.3, which is the average value of the Terlevich et al. (20011 
sampl^^ These transformations ensure that, when comparing all these samples with the LCBG models at 
z = 0, colors and magnitudes are free from the contribution of emission lines and the subsequent comparisons 
are between purely stellar properties only. 

It can be seen in figure |5] that the first set of models ends its evolution much closer to the bulk of the 
Sph systems than the second set of model does. In fact, should the metallicity of LCBGs be on the lower 
end of the range allowed by their line ratios, the simple passive path assumed here could not transform 
them into Sph galaxies. To illustrate this point, a third, 10 Gyr old stellar population was added to models 
^2. The mass of this population was arbitrarly set to be 50% of the dynamical mass. Two different metal 
contents were tried for this primeval population, Z ~ 0.0004 and Z ~ 0.0001 since higher metallicities are 
not justifiable in the framework of these metal-poor models. This population turned out to be unable to 
change the final B — R color to a great extent, affecting only in the scale of 0.05 magnitudes. This possibility 
was not explored further, then. It is concluded that, if these intermediate redshift LCBGs are to become Sph 
galaxies following a simple passive evolution, the first set of models is likely to be a better approximation 
to their stellar inventory. In order to reach to more robust conclusions on the metallicity of LCBGs and 
their descendants, it would be needed to compare the metallicity of LCBGs with that of Sph galaxies. 
However, the metal content of Sph galaxies is very difficult to determine spectroscopically because of their 
low surface brightnesses, and only broad ranges are reported in the literature, mainly based in photometric 



'^^Both studies dealt with the Coma Cluster, which is located at lOOMpc. 
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measurements (see, e.g Rakos et al. ( 2003[ ); Jerjen et al. (2004l). In the former work, the [Fe/H] metaUicity 



of dE galaxies is found to he in the 1/100 to 1.0 solar interval, whilst in the latter work the iron content of 
dE galaxies is found to range from 1/25 to 1/3 of the solar value. These broad ranges prevent us from using 
these values to discriminate between the presented models. If solar proportions are assumed, the previous 
iron content ranges roughly translate into 12 + logO/H ranging from 7.0 to 8.5. This range is unfortunately 
in agreement with both the metallicities of models #1 and #2. 

It is also seen in figure |5] that, even though BCDGs are less luminous than Sph systems or H II galaxies, 
their B — R colors are similar although somewhat bluer on average. BCDGs are as luminous as the brightest 
globular clusters, and their stellar populations are also vey old, like those of Sph systems. If they are the 
descendants of a hypothetical class of low-luminosity LCBG-like galaxies in the burst phase, the blue absolute 
magnitude of these objects should be around -16.0. Such objects would be very hard to observe and they 
would very similar to dwarf irregular systems like NGC1569 or NGC1705. 

Figure [5] shows that the underlying stellar populations of the local H II galaxies are compatible with the 
evolution of LCBGs in the sense that their B — R colors were the same at some point in the past around 
the time when the observed starbursts took place. It is then likely that the underlying stellar populations 
found in the presented LCBGs are much younger than the underlying stellar generations found in local HII 
galaxies. The underlying stellar population of local H II galaxies is found to be older that 1 Gyr, in the work 



of Westera et al. ( 2004 1 . The older stellar generations derived in the simple models presented gere resemble 



a somewhat aged intermediate age population, again in the nomenclature given in Westera et al. (2004) 



In order to compare to a greater detail the stellar populations of LCBGs and H II galaxies, the evolution 
in the {V — R)-{R — I) plane of the models devised here is presented in figure |6] together with the location 
of the underlying stellar populations of the 15 bright, high equivalent width H II galaxies from Telles & 



Terlevich (19971. These latter objects are again represented by half- filled dots. The V — R colors are in a 
very narrow range around 0.3, the R — I color has a much greater range of variability, from to 1. This 
figure also shows that the V — R color of the underlying population found in local H II galaxies is indeed 
attainable by LCBGs at some point during their evolutioij^ The underlying stellar populations of local H II 
galaxies are then found to emit more power in the I band than intermediate redshift LCBGs of similar V — R 
colors. This might be caused by the fact that intermediate redshift LCBGs probably harbor lower relative 
numbers of red giant stars than the redder H II galaxies since they have had much less time to evolve since 
the onset of the red giant phase of the evolution of their dominant stellar populations. Such very luminous 
red giant stars, existing in the extensions and wings of the local HII galaxies, would be responsible for a 
sizeable fraction of the I band luminosity of local H II galaxies. 

These suggestions can be interpreted in the light of the picture for H II galaxies presented in |Westera et al.| 



(2004) in which H II galaxies are a mixture of three populations. The first of them is the youngest population 
responsible for the observed emission lines, the second population is an intermediate age generation, older 
than 50Myr but younger than a few hundred Myr and the third population is the oldest one, older than IGyr. 
The single underlying stellar population used here can be thought to represent a mixture of the second and 



third populations of the Westera et al. (2004 ) model. As such, it can only accurately reproduce the colors of 



the real mixture it represents in a limited wavelength interval. Figure [6] shows that the model fails in the I 



'^^In particular, according to models #1, the underlying populations of LCBGs presented such V — R colors 4.9 Gyr ago, 
when the observed starbursts took place. In the case of the metal poor models, this happens 4.5Gyr ago, just after the observed 
starbursts occured. The first set of models reaches the average R — I oi the local sample 2.3Gyr ago, while models #2 will 
attain such R — I color in the future. 
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band. It is then natural to think that the real stellar mixture in H II galaxies is relatively richer in very old, 
I emitting stars, namely red giants. This very old and important stellar component has already been shown 
to exist in local H II galaxies. In addition to this, the work presented in Raimann et al. (20001 showed a 



detailed statistical analysis of the stellar populations found in a large sample of emission- line galaxies. It was 
found that a significant fraction of the stellar mass built in these systems is at least 500Myr old. Present day 
descendants of LCBGs will then have very old stellar populations, similar to the oldest stellar populations 
described in Westera et al. (20041. Should these systems undergo another star formation episode they would 



be very similar to modern H II galaxies. More modelling and observational work is clearly needed in order 
to ascertain this. 

Figure [7] presents other aspects of the evolution of the constructed models. It shows the evolution of 
the stellar B — Fcolor as a function of look-back time. The symbols at t ~ represent four well known 
spheroidal systems, as indicated in the figure. The metal rich models end their evolution near the position 
occuppied by the dE/Sph sources. It can also be seen that the second set of models does not reach the very 
red colors of the presented Sph galaxies. Moreover, these metal poor models can only reach B — V= 0.7 in 
the distant future. Such colors are not in agreement with the colors of the handful of Sph galaxies presented 
for comparison. 

In figure [s] the evolution of the presented models in the \og{M / Lb)^{B~V) diagram is shown. The 
positions of the LMC and individual dwarf elliptical systems is shown with specific symbols for each object. 



The position of the LCBGs observed with the WF/PC-2 in Guzman et al (1998 1 is depicted using small 



squares with asymmetric error bars. The errors in M/ L arise from the fact that the masses used in Guzman 
et al ( I998|) were not stellar masses but dynamical ones. The presented error bars allow for a factor of two 
difference in both numbers. The errors in the _B — y color come from the contribution of the emission lines 
to the spectrum. This contribution has been assumed to be 0.05 magnitudes towards bluer colors for all 
objects. The B — I/color of the LMC has been corrected for the contribution of the line-emitting region and 
starburst 30-Dor and extinction by subtracting 0.15 to the color. 

It is again seen that models #1 end their evolution much closer to the Sph galaxies, although neither 
model falls neatly within the region defined by the Sph galaxies shown. The fact that the mass-to- light ratios 
of the first set of models are more similar to those of the Sph systems than those of models #2 allows us to 
think that the rest of the other colors of dE will be more similar to the colors predicted by models #1. It is 
also the case that the stellar populations of the presented models have a mass-to-light ratio lower to that of 
the LMC by the time their B — V color is similar. This suggests that the star formation has proceeded for a 
longer period of time in the LMC than in LCBGs, which is again consistent with LCBGs being observed at 



very high look-back times. The Guzman et al (1998) objects cluster around the model tracks. In fact, six 



objects lie directly on top of the models themselves. This is reassuring in the sense that the bulk of the stellar 



inventory in the Guzman et al ( 1998 1 sources is probably well represented by the simple models presented 



here. However, it is somewhat unfortunate that the Guzman et al (1998 1 points fall both above and below 



the model lines because this makes it impossible to further resolve the underlying stellar population. If, for 



instance, the Guzman et al (1998 1 points were all above the model tracks around the LMC, it would have 



been possible to argue that the real underlying stellar population should be older than the age predicted by 
the models. This fact, however unfortunate, was to be expected since the i? — V^color is not very informative 
in the age range studied. The use of redder colors would be needed in order to gain this drop of information 
on the underlying stellar population. 



It is also interesting to compare figure |8] with Guzman et al (1998 Fig 3. of). It is seen that the path 
of the models built in this work in the \og{M / L b)~{B — V) plane is very similar to the path of the Guzman 



B-R 



B-R 



Fig. 5. — Model evolution in the {B — R)-Mb plane. The green line represents the evolution of the SA57-7042 
models, the blue line represents SA57-5482, the red line shows the evolution of SA57-10601 and the black 
line is Hl-13088. The black dots represent the dwarf elliptical systems from 'Guzm an | (19941. Half filled 
dots are local H II galaxies from Telles & Terlevich (19971. Crosses and white squares represent the local 



BCDG samples from Doublier et al. ( 1997[ 19991; Caon et al. (2005) respectively. The left panel represents 
the metal-rich models whereas the right panel depict models #2. 
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Fig. 6. — Both panels represent the evolution in the V — R-R — I plane of the models used here to describe 
the four observed LCBGs. The left panel represents the first set of models, while the right one shows the 
metal-poor models. The local data are taken from the HII galaxy sample from ,Telles fc Terlevichj ( |1997[ ). 
The square star represents the average position of this H II galaxy sample. Other symbols as in figure [Sj 
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et al (19981 models, implying that the differences in the age of the underlying stellar generation and relative 
mass of the ionizing population have little impact on the _B — V^color and Mass-to-Light ratio, at least within 
these limits. Other observables will surely be more affected. 

Figure |9]is very similar to figure [s] It shows how the presented models evolve in the Mb-{U—B) plane. 



This time, only the position of the LMC and of the sample of spheroidal systems taken from Guzman ( 1994 1 
is included. The H II galaxies and BCDGs are not included because there are not U band data for these 
galaxies, and estimating it from V band data would be affected by large uncertainties. The contribution to 
the C/— Scolor from the 30-Doradus cluster and extinction is assumed to be 0.1. It is again seen that the 
final stage of the first set of models is much more similar to the observed properties of the Sph galaxies. 
Furthermore, all the models converge towards a point around U—B = 0.3, Mb = —17.0. This suggests that 
any Sph descendant of the distant LCBGs should be the very red and luminous. This latter conclusion is 
weaker in figure |5] However, since the model used is very simple, it is also likely that the aforementioned 
convergence towards this subset of the Sph galaxies merely reflects a weakness of the models. Indeed, should 
any major star forming episode happen in these objects, the final position of the evolved LCBGs would be 
different. It is also interesting to note that the LCBG models run very near the LMC at some point in the 
past, some time after the observed starbursts took place. A larger sample of LCBGs with more data would 
be needed to study the reality of the convergence effect mentioned above. 



Finally, figure 10 shows the evolution of the stellar populations predicted by our models in the (U—B)- 
{B — R) plane. It is again seen that the first set of models ends its evolution much closer to the Sph galaxies 
than the metal poor ones. Moreover, the models converge towards the bulk of the dwarf elliptical galaxies, 
and this conclusion also holds when only the bright Sph systems are considered, strengthening the idea that 
the Sph descendents of LCBGs are the more luminous ones. 



6. Summary and Conclusions. 

We have used spatially resolved STIS long-slit spectroscopy together with optical and NIR imaging to 
investigate the stellar populations found in the the line-emitting region of four LCBGs at intermediate z. 
This is one of the first times this has been done for distant, compact galaxies. The star bursting episodes 
are found to be less than 3.5 Myr old in most cases, with typical masses in the range from 5 x 10^-12 x 10^ 
solar masses. The mass contribution of the newly formed stars to the total mass of the stellar component 
goes from 0.1% to 4%, with an average value of 1%. The ages of the older stellar generations found in these 
systems arc in the 20 Myr to 1400 Myr, with an average value of 700 Myr. The NIR observations were 
used to check for consistency the optically derived stellar populations. It turns out that the agreement is 
generally very good, implying that the real stellar content is bound to be very similar to the two population 
mixture used here. The line-emitting region spans ~ 2 kpc in size and is found to be very concentrated in 
the observed galaxies, although it is not always in the center of the optical galaxies. The metal content of the 
ionized gas phase can not be readily determined with the data at hand, although the nitrogen measurements 
of very similar objects suggest slightly subsolar abundances. The NICMOS data allow to regard the metal 
rich models as the favored ones, too. 

Although the available data are very limited, we have shown evidence for a very important underlying 
stellar population even in the central region of these sources, where the starburst completely dominates 
the observed optical emission. This population is somewhat evolved and accounts for most of the mass of 
these systems, despite the high look-back times probed. It can be said that the presented star formation 
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Fig. 7. — Temporal evolution of the stellar B — Vcolor with look-back time. The left panel represents the first 
set of models, while the right one shows the metal-poor ones. Individual symbols represent the B — t^color of 
several well known local spheroidal systems. The black pentagon is NGC205, the white triangle is IC3393, 
the cross is NGC3605 and the white circle is NGC147. Other symbols as in figure [5] 




Fig. 8. — Model evolution in the B — V-{M/Lb) plane. As before, the left panel represents the first set of 
models, while the other panel represents the metal-poor models. The position of the galaxies with M/L 
ratios taken from Guzman et al (19981 is given as small squares, and their error bars are explained in the 
text. Individual symbols represent the position of several well known local spheroidal systems. The Large 
Magallanic Cloud is included, too. The white triangle is IC3393, the big cross is NGC3605, the white circle 
is NGC147 and the black hexagon is the LMC. Other symbols as in figure |5] 
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Fig. 9. — Model evolution in the {U — B)-Mb plane. The Large Magallanic Cloud is included as a black 



hexagon. Both open and closed circles represent the Sph galaxies from |Guzman (1994). Closed symbols 
represent objects brighter than AIb = —15.0. Open symbols are the dimmer objects. 




Fig. 10. — Evolution in the [U — B)-{B — R) plane of the stellar populations that best represent the observed 
LCBGs. Both open and closed circles represent the Sph galaxies from Guzman (1994). Closed symbols 
represent objects brighter than Mb = —15.0. Open symbols are the dimmer objects. Panel layout as in 
figure (5) 
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histories are similar to some extent to the star formation histories of HII galaxies described in Westera et al. 



(|2004) . Both HII galaxies and LCBGs are therefore age-composite systems. They both have an old stellar 
population, together with a young stellar generation. It has also been found that, should LCBGs finally 
evolve into dwarf ellipical galaxies, their metal content must be in the higher end of the interval allowed 
by their observed line ratios. Unfortunately, given the broad range of the metallicity determinations of Sph 
systems, this information can not be used to help rule out any one of the presented models. The NIR data 
point towards the higher metallicity ones, though In addition, it's been concluded that the Sph descendants 
of LCBGs can only be the more luminous and red of such systems. The comparison of the LCBG models 
with the underlying population of luminous and local HII galaxies indicates that the star formation histories 
of the observed LCBGs and local H II galaxies are similar. It has to be kept in mind, however, that H II 
galaxies probably have a larger relative number of red giant stars, as hinted by their enhanced power in the I 
band. This can be simply interpreted by the fact that they are, arguably, 5 Gyr older than the more distant 
LCBGs. The observed differences in the redder colors probably arise from differences in the relative numbers 
of red giant stars. This suggests that the recent star formation histories of LCBGs and H II galaxies are 
similar up to some point. From that point backwards, the star formation histories probably differ. It has 
also been pointed out the possible existence of dwarf irregular systems in the intermediate redshift range. 
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